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Abstract 

We present the one-loop corrections originating from Quantum Chromo-Dynamics 
(QCD) and Electro-Weak (EW) interactions of Supersymmetric (SUSY) origin 
within the Minimal Supersymmetric Standard Model (MSSM) to the single-top pro- 
cesses bq — )• tq' and qq' — )• tb. We illustrate their impact onto top quark observables 
accessible at the Large Hadron Collider (LHC) in the 't+jet' final state, such as total 
cross section, several differential distributions and left-right plus forward-backward 
asymmetries. We find that in many instances these effects can be observable for 
planned LHC energies and luminosities, quite large as well as rather sensitive to 
several MSSM parameters. 
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1 Introduction 



Top quark processes at the LHC will be an ideal laboratory where it will be possible 
to profile with great accuracy the heaviest particle of the Standard Model (SM) and 
eventually extract possible effects of New Physics (NP) Beyond the SM (BSM). The 
manifestation of NP at the LHC will in fact depend on the actual NP mass scale, Mnp- 
If the available energy is less than Mnp, then NP will manifests itself through virtual 
effects. Conversely, if the available energy is greater than Mnp, NP will appear through 
the production and decay of new real particle resonances. A mixture of both scenarios 
may of course also occur. 

If NP is identified with SUSY, then it is fair to say, based on the available literature, 
that much has been done in the second respect. Various methods to extract the presence 
of new resonant SUSY particles have in fact been devised over the years and applied 
to several possible model realisations of a SUSY theory, particularly for the case of the 
simplest SUSY scenario, the MSSM. In constrast, much less has been done so far in the 
first respect. Primarily because the LHC is perceived not to be a precision machine, due 
to both the undefined partonic energy and the intrinsically large hadronic background, 
so that virtual effects (of order a.^ and a or indeed smaller) are expected to be not easily 
discernible in the data. However, there are certain processes where the experimental 
precision is expected to eventually become comparable to the size of the virtual effects of 
NP. 

Among these processes, one can certainly list top quark production, both in double- 
and single-top mode. In fact, with typical hadro-production cross sections at the LHC of 
order 800 pb and 300 pb respectively (at 14 TeV) and collider luminosities that can reach 
the 300 fb~^ level, several hundred miUion top events will altogether be produced during 
the lifetime of the CERN machine. Clearly, on the one hand, this renders the statisti- 
cal error applicable on the experimental side to typical top quark observables negligible. 
Furthermore, on the other hand, the main systematics affecting the latter (both theoreti- 
cal, coming from the Parton Distribution Function (PDF) dependence, and experimental, 
coming from the jet energy scale uncertainty) will be understood at the percent level after 
the first hadronic data samples will have been collected and studied in detail. 

Therefore, both double- and single-top production processes at the LHC lend them- 
selves to being precision physics laboratories, where, in particular, virtual effects of SUSY 
(when Mnp = Msusy ^ iTit, where rrit is the top quark mass) can possibly be extracted, 
in the ultimate attempt to understand the underlying dynamics of SUSY breaking. We 
intend to illustrate here the potential afforded in this respect by two single-top production 
processes, namely the t-channel: 

bq tq' (1) 

and the s-channel: 

qq' tb (2) 
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(plus, of course, their charge conjugated channels), where g*^'^ represents a light quark 
(m, d, s or c) 

The plan of the paper is as follows. The next section describes what is available in 
current literature concerning single-top processes at the LHC in terms of their higher 
order corrections within the MSSM. Sect. [3] illustrates how we performed the calculation 
of the one-loop SUSY QCD and EW corrections to the two single-top channels ([I])-(l2]). 
The following section presents our numerical results. We conclude in Sect. |5l 

2 Single-top Processes at the LHC 

Due to the their relevance for LHC physics, single- as well as double-top processes have 
extensively been studied in previous years, with a twofold aim: on the one hand, in order 
to have a precise SM prediction, at least at the complete NLO (and possibly next-to-NLO 
(NNLO) for the QCD corrections); on the other hand, in order to investigate possible 
deviations from SM predictions due to the presence of NP. These two goals are obviously 
related, the former being a precondition for the second: in order to extract meaningful 
information from NP processes entering at one-loop level top quark hadro-pro duct ion at 
the LHC, clearly, all similar SM effects should be well under control. 

One-loop SM corrections to single- and double-top production in hadron-hadron col- 
lisions have been known for some time. These include the SM QCD pLj and SM EW |2] 
ones to gg, qq — )■ ti (double-top production) and the corresponding ones [Hill] for bq — )■ tq', 
qq' — > th and hg — )■ tW~ (single-top production). Both SM QCD and EW corrections can 
be large, although in complementary energy regimes: ^/s ~ mt and ^/s ^ m^, respec- 
tively. However, they all are rather stable against variations of the factorisation and/or 
renormalisation scales (these in turn quantifying the systematic uncertanties related to 
the unknown two- loop, or even higher order, corrections), so that one can conceivably 
attempt to investigate virtual effects of some NP, particularly of SUSY origin, induced 
onto single- and double-top processes. To stay with the MSSM, QCD and EW corrections 
to double-top production have been quantified in [5] and [6], respectively, whilst only part 
of these are known to date for the case of single-top reactions [7]. 

Our paper is thus the first one where both MSSM QCD and EW corrections to single- 
top processes inducing a 't plus jet' signature in the final state are computed (i.e., hq — )■ tq' 
and qq' — )■ th). It is in fact from this perspective that we decided to postpone to another 
publication the presentation of similar corrections to the third single-top channel (i.e., 
hg — )■ tW~), in the sense that the latter originates an altogether different signature in the 
final state, 't plus lepton and missing (transverse) energy' or 't plus two jets', depending on 
whether the boson produced in association with the top quark decays leptonically or 
hadronically. In fact, both of the latter require different triggers and selection procedures 
with respect to the former, so that, from a phenomenological point of view, they deserve 
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a separate treatment. 

3 Calculation 

The values of the SM input parameters considered for the numerical evaluation of the 
one-loop corrections can be found in Tab. [TJ 



Coupling constants 


a{Mz) = 1/127.93400652 a,(Mz) = 0.1176 


Gauge boson masses 


Mw = 80.424 GeV Mz = 91.1876 GeV 


Heavy Quark Masses 


mt = 170.9 GeV mf^irrib) = 4.2 GeV 



Table 1: Numerical values of SM inputs. 



We used the CTEQ6 2006 PDF set [8J, in particular the CTEQ6L fit, which is LO in 
QCD, choosing renormalisation and factorisation scales both equal to the top mass, m^. 

As intimated in the introduction, due to the relatively large single-top production 
cross sections (in particular for the t-channel process) and the consequent large available 
statistics, we reckon that it will be possible at the LHC to perform "precision physics" of 
top samples produced singly. It is therefore a worthy task trying to understand system- 
atically if sizable virtual effects could modify the SM predictions, providing an indirect 
but encouraging hint of NP. As a first step we performed an adaptive scan over both the 
MSSM and minimal- S Up er GRAvity (mSUGRA) parameter spaces in order to investigate 
whether there exist regions with large one-loop corrections. We have done so limitedly to 
the inclusive production cross sections in either channel. 

While in the MSSM case we have input the relevant parameters directly at the EW 
scale, in the mSUGRA scenario the low energy spectra have been obtained evolving the 
input parameters from the Grand Unification Theory (GUT) scale down to the EW scale 
through the code SUSPECT [9]. 

We have developed a dedicated C++ code to compute the one-loop corrections and, 
as an internal check, we have tested the cancellation of Ultra- Violet (UV) divergences 
appearing in the loop integrations. To quantify the one-loop corrections, we will mainly 
focus on ii'-factors, defined in general as the ratio NLO/LO of a given observable, with all 
relevant quantities (QCD and EW coupling constants, etc.) evaluated at the given order 
consistently. The PDFs, on the other hand, have been kept at LO since there NLO effects 
are of SM origin, whereas we are looking here at purely SUSY corrections, therefore using 
NLO QCD PDF sets would introduce spurious SM effects from higher orders. 

Guided by the inclusive results, we have then looked at differential distributions for 
some 'benchmark points'. In particular, alongside the typical kinematical observables 
(invariant masses, transverse momenta, (pseudo)rapidities), due to the possibility of mea- 
suring spin-related observables in top samples (the top (anti) quark decays before hadro- 
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nising in fact, thus efficiently transmitting its polarisation state onto the ensuing decay 
products), we have also focused our attention on the following observables: 



• the Left-Right Asymmetry (Alr); 

• the Forward-Backward Asymmetry (A-pb), only for the s-channel. 
Alr is defined as usual by the following ratio: 

-^LR = ■ • [O) 

Cpp^tL+X + CTpp^tn+X 

ApB deserves further details, since the definition of "Forward" and "Backward" re- 
gions in a pp collider is not straightforward. In the s-channel process the initial state 
can essentially be ud or cs. In the first case m is a valence quark, and thus the momen- 
tum fraction it carries is much bigger than that of the d sea-quark. It is then possible 
to reconstruct here the direction of the incoming valence quark just by identifying the 
direction of the boosted top in the final state. This is a statistical process, since s-channel 
single-top production can be initiated by two sea-quarks too, such as in the cs case, and 
the weight of the various events is given by the relevant PDF (hereafter denoted as P). 
The definition of A-pB is then given by the following relation [10] : 



/ dXidX2 'Eq=u,d,cA^1^^^' f^)Pg{X2, + Pg{Xi, fi)Pg{X2, /i) ^ ^^^^ 

dxidx2(T{pp ^ tb + X) (5) 
{Pu{xi,fx)Pa{x2, /i) - Pd{xi,n)Pu{x2, /i))sign(xi - X2)] (6) 

where fi is the factorisation/renormalisation scale. 

We also have computed differential distributions of these quantities over the invariant 
mass of the final state, the transverse momentum of the top as well as the rapidity of 
both the top and the light quark in the final state. 

Coming back to the scans, in the general MSSM scenario, we have explored over all 
soft parameters involved in the generation of the spectra, while in the mSUGRA scenario 
we have scanned over the four standard mSUGRA input parameters. The scans have 
been performed interfacing the tool adScan[Tl] to our code. 

The main result of this preliminary analysis, both for MSSM and mSUGRA, is that 
the relative effect of the SUSY one-loop corrections seems to be rather small over the 
whole range of the independent parameters. Moreover, the effects show a very smooth 
and mild dependence on the parameters. A meaningful sample of our results is shown 
in Fig. [TJ Each point in the diagrams corresponds to a full computation of NLO effects 
for a point in the parameter space, plotted versus a low energy parameter, such as a 
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mass or mixing angle evaluated at the EW scale. We have chosen to show just three 
of them, corresponding to the parameters which are expected to affect most heavily the 
NLO corrections: rrig (the gluino mass), tan /3 (the ratio of the two doublet Higgs vacuum 
expectation values) and m^^ (the lightest stop mass). 

From the tan /3 plot it is easy to see that the density of the points is larger for very small 
values of the corrections and rapidly decreases for larger values of the effect. Furthermore, 
and possibly contrary to naive expectations, the density of the points is almost tan /3 
independent. In essence, for the largest part of the parameters space, the corrections 
are very small, below the 2% level, while it is unlikely to find parameters configurations 
leading to corrections as high as 5%, and this result holds independently of tan/3 (i.e., 
there are no visible regions in the tan /3 range where one-loop effects have a maximum). 
Trivially, the plots of the corrections as a function of the gluino or stop masses simply 
state that it is more likely to find "larger" (6% or so at the most) corrections in parameters 
regions where these masses are small. 

Despite the results obtained from the scan are never very large, it should be noted 
that these refer to inclusive corrections only. However, while scanning over the MSSM 
and mSUGRA parameter spaces, we have also tested for the size of SUSY QCD and EW 
corrections to differential distributions (including asymmetries) and found them to be 
rather large, albeit in limited (and at times disfavoured) regions of phase space, thereby 
explaining their smallness in the inclusive results. In order to illustrate the typical pat- 
tern emerging at differential level, we have defined a couple of representative benchmark 
mSUGRA points, whose parameters are given in Tab. [2l 



mSUGRA scenario 


mo 


mi/2 


Ao 


tan P 


sign /i 


LS2 


300 


150 


-500 


50 


+ 


SPSla 


100 


250 


-100 


10 


+ 



Table 2: Input parameters for the mSUGRA benchmark points (all values with mass 
dimension are in GeV). 

The SPSla point [12] is the typical mSUGRA "standard candle" while the LS2 point 
has been introduced in [13] and it is characterised by an interestingly light spectrum (yet 
compliant with current experimental and theoretical constraints), which - if realised in 
Nature - could prompt for easily detectable signals at the LHC. More details about the 
two spectra can be found in Fig. |2l 

4 Results 

It has recently been scheduled a long LHC run at 7 TeV, which will probably last two 
years and will hopefully lead directly to the final expected energy of 14 TeV. In view 
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of this recent development, we present our results for both such planned energies. Our 
numerical results will be illustrated in turn in the following two sub-sections. Notice that 
we will not show all the differential distributions analysed for the two benchmark points, 
rather we will illustrate a meaningful sample of these (and only in the 14 TeV case) and 
focus on the one-loop corrections they receive. We will also show the contribution of the 
SUSY EW and QCD parts separately, in some instances at least, to better understand 
the features of the underlying dynamics: like, e.g., dominant contributions, asymptotic 
behaviours and threshold effects. 

4.1 LHC at 7 TeV 

The total cross section and asymmetries for t- and s-channel are given in Tabs. |3] and H] 
respectively. 





a(pb) 


L-R Asymmetry (1-5) 


Born 


35.44 


(5=4.7327x10^^ 


LS2 
SPSl 


35.29 (ir=0.996) 
35.33 {K=0.997) 


5=4.7515x10^^ (l-is:=l. 87x10"^) 
5=4.7353x10-^ (l-fs:=2.6xl0-i°) 



Table 3: Numerical results for t-channel production at 7 TeV. 





(T(pb) 


L-R Asymmetry 


F-B Asymmetry 


Born 


2.061 


0.6777 


0.598049 


LS2 
SPSl 


2.086 (Js:=1.0121) 
2.062 (ir=1.0005) 


0.6796 (ir=1.0028) 
0.6783 (ir=1.0009) 


0.598018 (fs:=0.999948) 
0.598046 (fs:=0.999995) 



Table 4: Numerical results for s-channel production at 7 TeV. 



As it is possible to see, in both the two channels, the NLO corrections to the total 
cross sections are opposite in sign: in the t-channel case one-loop corrections reduce the 
cross section, while in the s-channel there is an enhancement, but the corrections are 1% 
or less in both the benchmarks considered (in line with the results from the mSUGRA 
scans) . The main difference between the two channels is of course the value of the cross 
section, a. Focusing on t-channel, which has a greater a, we have found that, with an 
integrated luminosity of 1 fb~^ the LHC will generate ~ 35000 events per year, therefore 
the detection of any deviation from the SM prediction due to the SUSY EW and QCD 
one- loop corrections, at least in the early stages of the experiment, will be indeed very 
challenging. Also Alr shows different behaviours in the two channels but, unfortunately, 
they are much greater in the s-channel, which, as we have just verified, has NLO correc- 
tions well beyond observability at 7 TeV. The same observability considerations can be 
applied to A-pB in the s-channel. 
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4.2 LHC at 14 TeV 
4.2.1 The t-channel case 

The total cross section and asymmetries for this production mechanism are given in Tab. |5l 





(T(pb) 


L-R Asymmetry (1-5) 


Born 


122.5 


5=4.3491 xlO"^ 


LS2 
SPSl 


122.0 {K=0.996) 

122.1 (i^=0.997) 


5=4.3668x10"^ (l-i^=l. 77x10-9) 
5=4.3516x10-7 (l-ir=2. 5x10-1°) 



Table 5: Numerical results for t-channel production at 14 TeV. 



In this case the total cross section is significantly bigger than in the 7 TeV case, 
and therefore we expect that it will be possible to extract much more information from 
experimental data. The i^-factor for the integrated cross section is admittedly small, less 
than 1% for both benchmarks but, with the expected luminosity of 10 fb-^, the production 
of single-top through this t-channel process will count ~1M events and therefore it should 
be possible to detect a discrepancy in the expected number of SM events that could be 
interpreted as a hint of NP and thus boost the search of a SUSY resonance. As far as 
the Alr is concerned, the top quark will emerge almost completely left-polarised, yet the 
detection of a one-loop induced discrepancy in the asymmetry will be again unlikely, the 
X-factor being close to one part in a billion or less. 

The analysis of differential observables is now mandatory to better understand the 
origin of the corrections to the integrated cross section and Left-Right asymmetry. In 
Fig. [3] we have collected our results for the differential distributions of the cross section. 
It can be noticed that the corrections are always bigger in the LS2 scenario than in the 
SPSla case and that in the transverse momentum distribution they can reach the —10% 
limit for high values of the top pt (transverse momentum). However, the differential cross 
section is quite small in this limit, of the order of 0.01 pb/TeV, and therefore observation 
of large one-loop corrections for events with high top transverse momentum is probably 
possible though quite difficult. It is also possible to see that the one-loop corrections in 
the other differential distributions are always quite small, below the 1% or so limit. 

In Fig. m the one-loop corrections to the differential distributions for the Left-Right 
asymmetry are shown. Again, the corrections are generally bigger in the LS2 scenario 
than in the SPSla case and can reach significant values for high top pt, of —8% or so. 
The Pt distribution has another interesting feature: the asymmetry is bigger in the high 
top Pt region, where, however, one-loop corrections are much less than 1%. 

4.2.2 The s-channel case 

The total cross section and asymmetries for this production mechanism are given in Tab. IB 
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cT(pb) 


L-R Asymmetry 


F-B Asymmetry 


Bom 


5.138 


0.6861 


0.53665 


LS2 
SPSl 


5.206 (ir=1.0132) 
5.145 (^=1.0014) 


0.6883 (ir=1.0032) 
0.6869 (ir=1.0012) 


0.53695 (fs:=l. 00056) 
0.53676 (fs:=l. 00020) 



Table 6: Numerical results for s-channel production at 14 TeV. 



Single-top events from s-channel at the LHC will be significantly more than at the 
Tevatron, and this will allow a much more precise analysis of the process and of deviations 
to related observables due to new physics. The inclusive one-loop corrections we have 
obtained are at most of ~1% in the LS2 scenario, but maybe enough to allow a small 
observable difference in the expected number of events. In constrast, the Left-Right 
asymmetry receives NLO corrections of the order of 0.1% in both the scenarios considered 
(here, unlike the general case, the effect of pure SUSY corrections is almost completely 
negligible). The same result holds for the one- loop corrections to the Forward- Backward 
asymmetry: the ii'-factor is very small and the effect practically undetectable with the 
s-channel cross section being of some pb. 

The differential distributions for the cross section are shown in Fig. [51 In the s-channel 
case there are significant differences between the two scenarios considered: the corrections 
for LS2 are higher than those for SPSla in the regions where the contributions to the total 
cross section are bigger, and this is the reason for the big difference in the i^-factors for 
the integrated cross section in the two scenarios. 

The same feature can be observed for the Left-Right asymmetry distributions. Fig. |6l 
where, however, the enhancement of the one-loop corrections from the LS2 scenario is 
milder and therefore the differences in the integrated asymmetry are smaller. 

The corrections to differental quantities in the Forward-Backward asymmetry are qual- 
itatively and quantitatively very similar to those for the Left-Right asymmetry, as it is 
possible to see in Fig. U\ thus the same comments as above apply in this case. 

Considering all the differential distributions for the s-channel process, we can notice 
that they do not differ too much from each other between the two benchmarks and, as 
a general feature, they reach their maximum value in the invariant mass distribution for 
Minv ^ 700 GeV, in the transverse momentum distribution for p-r > 300 GeV and in the 
top rapidity for yt > 1.5. 

To better understand the origin of the corrections and why there are peaks for cer- 
tain values of M[^^, pt and yt we have computed (here, limitedly to the s-channel) the 
SUSY EW and QCD corrections for the cross section separately and plotted the related 
differential quantities in Figs. El and [9] for the two benchmark scenarios. It is possible 
to see that the corrections are dominated by the SUSY QCD contribution, while the 
EW counterpart provides small contributions to the total corrections. In both bench- 
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mark scenarios, however, the EW contribution shows pecuhar features in the invariant 
mass distribution, given by peaks and troughs in the correction, and a closer inspection of 
these threshold effects reveals that they are situated in correspondence of . Therefore 
charginos could play an interesting role in the determination of SUSY effects in s-channel 
single top production processes. Going back to the dominant QCD contribution, the only 
parameter which enter the QCD correction alone is the gluino mass {rrig = 607 GeV in 
SPSla and rrig = 392 GeV in LS2), thus we can argue that the one-loop corrections are 
sensitive mostly to this parameter. Smaller gluino masses shift the SUSY QCD peak in 
the corrections towards regions where the differential distributions have higher values and 
therefore the integrated quantities will be affected more by NLO corrections. This means 
that observing higher pure SUSY one-loop corrections would point towards scenarios with 
lighter gluino and viceversa. 

To conclude this section, we are aware that recent measurements have slightly changed 
the value of the top mass|14j. Nevertheless, we have verified that any change (within a 
reasonable range) in this value has practically no effect on one-loop corrections, as it is 
possible to see in Fig. [10 

5 Conclusions 

In summary, QCD and EW corrections through one-loop level originating from SUSY 
and affecting the cross section for 't -|- jet' final states at the LHC are always small at 
inclusive level, at both energies of 7 and 14 TeV. However, they appear sizable (up to 
ten percent) and detectable (after (9(100) fb~^ of luminosity) at differential level, if the 
collider energy is 14 TeV, while for 7 TeV they will presumably not be accessible. Such 
large corrections also appear in the Left-Right and Forward-Backward (when defineable) 
asymmetries, again, at the differential but not inclusive level. The bulk of the higher order 
effects comes from QCD, with a much smaller, yet very sensitive to SUSY parameters, 
EW component. The corrections discussed here are typically larger for the s-channel mode 
than for the t-channel one, which renders them more difficult to cLCCGSS cLS the former is 
subleading with respect to the latter at the LHC. Altogether though, our results point to 
the relevance of SUSY one-loop effects onto 't + jet' final states for the LHC running at 
design energy and luminosity. In presence of high statistic samples, some typical SUSY 
parameter dependencies could possibly be extracted, in either mSUGRA or the MSSM. 
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Figure 1: Scan over three relevant SUSY parameters (tan/3, rrig and mj) in the context of 
the MSSM (upper panels) and mSUGRA (lower panels) to search for maxima of one-loop 
effects. 
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SPSla spectrum 




LS2 spectrum 




Figure 2: The SPSla and LS2 mass spectra. 
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Figure 3: Differential distributions of the cross section in t-channel at 14 TeV. In this 
and in the following Figs., the small panels focus on the range where the one-loop ef- 
fects are higher, and they are shown only when the corrections from LS2 and SPSla are 
distinguishable. 
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Figure 4: Differential distributions of the Left-Right asymmetry in t-channel at 14 TeV. 
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Figure 5: Differential distributions of the cross section in s-channel at 14 TeV. 
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Figure 6: Differential distributions of the Left-Right asymmetry in s-channel at 14 TeV. 
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Apg - Bottom Rapidity Distribution 
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Figure 7: Differential distributions of the Forward-Backward asymmetry in s-channel at 
14 TeV. 
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Figure 8: SUSY EW and QCD corrections to differential distributions of the cross section 
in s-channel at 14 TeV for LS2. 
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Figure 9: SUSY EW and QCD corrections to differential distributions of the cross section 
in s-channel at 14 TeV for SPSla. 
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Figure 10: K-factor dependence on the top mass in s-channel at 14 TeV for SPSla. 
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